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■ Abstract 

Qh! A detailed analysis is presented of the diffractive deep-inelastic scattering process ep — > 

eXY , where Y is a proton or a low mass proton excitation carrying a fraction 1 — x p > 0.95 
of the incident proton longitudinal momentum and the squared four-momentum transfer at 
the proton vertex satisfies \t\ < 1 GeV 2 . Using data taken by the HI experiment, the 
cross section is measured for photon virtualities in the range 3.5 < Q 2 < 1600 GeV 2 , 
triple differentially in x F , Q 2 and (3 = x/x p , where x is the Bjorken scaling variable. 
At low x p , the data are consistent with a factorisable x F dependence, which can be 
described by the exchange of an effective pomeron trajectory with intercept a p (0) = 
1.118 ± 0.008 (exp.) [to'oio ( m °del). Diffractive parton distribution functions and their 
uncertainties are determined from a next-to-leading order DGLAP QCD analysis of the Q 2 
and f3 dependences of the cross section. The resulting gluon distribution carries an inte- 
grated fraction of around 70% of the exchanged momentum in the Q 2 range studied. Total 
and differential cross sections are also measured for the diffractive charged current process 
e + p — ► v e XY and are found to be well described by predictions based on the diffractive 
parton distributions. The ratio of the diffractive to the inclusive neutral current ep cross 
sections is studied. Over most of the kinematic range, this ratio shows no significant de- 
pendence on Q 2 at fixed x p and x or on x at fixed Q 2 and (5. 
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1 Introduction 



Quantum Chromodynamics (QCD) is well established as the gauge field theory of the strong 
interaction. However, it is only able to provide reliable predictions for scattering processes if 
they involve short distance partonic interactions, where perturbative methods may be applied. 
In contrast, hadronic scattering cross sections are dominated by soft interactions, to which per- 
turbation theory is not applicable. In a large fraction of these soft interactions, often termed 
'diffractive', one or both of the interacting hadrons remains intact. Such processes are com- 
monly discussed in terms of exchanges with net vacuum quantum numbers, though the exact 
nature of these exchanges is not well known. 

The observation of high transverse momentum jet production in diffractive pp scattering [1] 
introduced the possibility of understanding the diffractive exchange in terms of partons. The 
presence of processes of the type ep — > eXp (figure 1) in deep-inelastic scattering (DIS) at 
low Bjorken-x at the HERA collider [2] offers a uniquely well controlled environment in which 
to study the QCD properties and structure of diffraction. Several measurements of the semi- 
inclusive cross section for this 'diffractive DIS' process have been made by the HI [3-6] and 
ZEUS [7-11] collaborations. 



(a) 
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Figure 1: Schematic illustration of the neutral current diffractive DIS process ep — > eXp, 
proceeding via virtual photon exchange. The dotted lines in (a) and (b) show the points at 
which the diagram can be divided under the assumptions of QCD hard scattering collinear 
factorisation and proton vertex factorisation, respectively. The kinematic variables defined in 
section 2 are also indicated in (a). 

The detailed explanation of hard diffraction has become a major challenge in the develop- 
ment of our understanding of the strong interaction at high energies and low x values [12]. A 
wide variety of models has been put forward to interpret the dynamics of diffractive DIS as 
well as its relationships to inclusive DIS and to diffractive hadron-hadron scattering [13-19]. A 
general theoretical framework is provided by the proof [20] of a hard scattering QCD collinear 
factorisation theorem [21-23] for semi-inclusive DIS cross sections such as that for ep — > eXp. 
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As illustrated in figure la, this theorem implies that the concept of 'diffractive parton distri- 
bution functions' (DPDFs) [22,24] may be introduced, representing conditional proton parton 
probability distributions under the constraint of a leading final state proton with a particular 
four-momentum. Empirically, a further factorisation has been found to apply to good approxi- 
mation, whereby the variables which describe the proton vertex factorise from those describing 
the hard interaction [5,6], as illustrated in figure lb. According to this 'proton vertex' factori- 
sation, the shape of the DPDFs is independent of the four-momentum of the final state proton. 
The dependence of the DPDF normalisation on the proton four-vector can be parameterised 
conveniently using Regge asymptotics, which amounts to a description of diffraction in terms 
of the exchange of a factorisable 'pomeron' (IP) [25] with universal parton densities [26]. 

Several authors have analysed diffractive DIS data to extract DPDFs [5, 10, 16, 27-32]. 
In [5, 6], HI measurements are studied, with the conclusion that the data are compatible with 
proton vertex factorisation at low fractional proton energy losses, x F . At larger x p , a sep- 
arately factorisable sub-leading exchange (R), with a different x F dependence and partonic 
composition, is present. The DPDFs extracted in [5] are dominated by the gluon distribution. 
Further tests of the factorisation properties of diffractive DIS have been made by comparing 
predictions using these DPDFs with hadronic final state observables such as diffractive jet [33] 
and heavy quark [34] cross sections. These tests have shown a remarkable internal consis- 
tency within the HERA DIS data. In contrast, the DPDFs extracted in DIS are not expected to 
be directly applicable to hadron-hadron scattering [20-22,35]. Indeed diffractive factorisation 
breaks down spectacularly when the DPDFs from [5] are applied to diffractive pp interactions at 
the Tevatron [36]. However, with the introduction of an additional 'rapidity gap survival proba- 
bility' factor to account for secondary interactions between the beam remnants [37], the HERA 
DPDFs remain an essential ingredient in the phenomenology of diffraction at the Tevatron and 
the LHC [38]. 

In this paper, a new measurement 1 of the diffractive neutral current DIS cross section is pre- 
sented. This is based upon HI data for which there is an absence of hadronic activity in a large 
rapidity region extending close to the outgoing proton beam direction. It is thus complementary 
to measurements such as [6], in which the leading proton is detected and measured. The proton 
vertex factorisation property is tested and the dependence of the diffractive cross section on x p 
is expressed in terms of an effective pomeron intercept a F (0). The dependence on x and Q 2 
is interpreted through a QCD analysis using the DGLAP [39] evolution equations at next-to- 
leading order (NLO) [40], from which new DPDFs are determined. The kinematic range of 
validity of the DPDFs is tested systematically and, for the first time, an assessment is made of 
the experimental and theoretical uncertainties. A first measurement is also made of the diffrac- 
tive charged current cross section, which is compared with a prediction based on the DPDFs 
extracted from the neutral current data. The ratio of diffractive to inclusive cross sections and 
its kinematic dependences are also investigated. Section 2 introduces the formalism adopted 
for the paper. Section 3 describes the cross section measurement. The results are presented in 
sections 4-7 and are followed by a summary in section 8. 

'The results presented here are in agreement with the previous HI measurement using similar techniques [5] 
throughout most of the measured kinematic range. The new measurement and DPDFs supersede the old due to the 
improvements in precision, kinematic coverage, theoretical modelling and detector understanding. 
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2 Diffractive DIS Kinematic Variables and Observables 



The data studied here are subsets of inclusive HI neutral and charged current DIS samples, 
arising from the processes e + p — > e + X' and e + p — > P e X', where the positron (with four- 
momentum k) couples to an electroweak gauge boson (q), which interacts with the proton (P). 
The usual DIS kinematic variables are defined as 

q=-q, * = 2p- q ' y = p— k ' w 

where Q 2 is the boson virtuality, x is the longitudinal momentum fraction of the proton carried 
by the struck quark and y measures the inelasticity of the process. The squared invariant masses 
of the positron-proton and gauge boson-proton systems are s = (k + P) 2 and W 2 — (q + P) 2 , 
respectively. 

The hadronic final state of any DIS event may be broken down into two systems X and 
Y, separated by the largest gap in the rapidity distribution of the hadrons relative to an axis 
defined by the exchanged boson and the proton in their centre of mass frame [5]. If the masses 
M x and M Y of these two systems are small compared with the mass W of the full hadronic 
final state, the two systems are expected to be separated by a large rapidity gap and a colourless 
exchange of well defined four-momentum may be considered to have taken place between them. 
As illustrated in figure la, the longitudinal momentum fractions, x p of the colourless exchange 
with respect to the incoming proton, and f3 of the struck quark with respect to the colourless 
exchange, are then defined by 

q-P 2q-(P-p Y ) 

Here, p Y is the four-momentum of the Y system and f3x p = x. The squared four- momentum 
transferred at the proton vertex is 

t=(P- PY ) 2 . (3) 

The rapidity gap selection (section 3.4) implies that the data analysed in this paper are domi- 
nated by the case where Y is a lone proton and \t\is relatively small. However, since the system 
Y is not detected directly, a small admixture of proton excitations and other systems such as 
leading neutrons is also present (see section 3.5). 

The neutral current data are presented in the form of a 'diffractive reduced cross section' 
cr^ 3 \ integrated over the ranges of t and M Y specified in section 3.5 and related to the differ- 
ential cross section measured experimentally by 

" Y + .a^\x F ,x,Q 2 ), (4) 



dx F dx dQ 2 xQ 4 

where Y + = 1 + (1 — y) 2 . Similarly to inclusive DIS [41], the reduced e + p cross section 
depends on the diffractive structure functions F 2 D( ^ and in the one-photon exchange 

approximation according to 

a D(3) = F D(3) _ |^(3). (5) 



For y not too close to unity, a r = F 2 l ' ; holds to very good approximation. In previous 
measurements of inclusive diffractive DIS at HERA, the data were presented in terms of F 2 D ® 
instead of <Jr^. 

Due to the smaller available data sample, the charged current measurements must be inte- 
grated over some or all of the kinematic variables. They are presented as a total cross section 
and single differentially in either x p , (3 or Q 2 . 



3 Experimental Procedure 
3.1 HI Apparatus 

The HI coordinate system is defined such that the origin is at the nominal ep interaction point 
and the polar angle 9 = corresponds to the direction of the outgoing proton beam. The region 
9 < 90°, which has positive pseudorapidity rj = — In tan 9/2, is referred to as the 'forward' 
hemisphere. 

The interaction region is surrounded by the central tracking system, which consists of a 
silicon vertex detector, drift chambers and multi-wire proportional chambers, all located within 
a solenoidal magnetic field of 1.15 T. The trajectories of charged particles are measured in the 
range— 1.7 < rj < 1.7 with a transverse momentum resolution of ct(pt)/pt — 0.005 px /GeV© 
0.015. The Forward Tracking Detector and the Backward Drift Chamber (BDC) provide track 
segments of charged particles at smaller and larger 9 than the central tracker, respectively. 

A highly segmented Liquid Argon (LAr) sampling calorimeter, covering the range — 1.5 < 
rj < 3.4, surrounds the tracking chambers and consists of electromagnetic and hadronic sec- 
tions. The energy resolution is a(E)/E ~ 11% /^/E/ GeV for electrons and a(E)/E ~ 
50%/ \[E~I GeV for hadrons, as obtained from test beam measurements [42]. The backward 
direction (—4.0 < i] < —1.4) is covered by a lead / scintillating fibre calorimeter (SpaCal), 
which also has both electromagnetic and hadronic sections. The energy resolution for electrons 
is cr(E)/E ~ 7%/y/E/GeV. 

In addition to the 'central detectors' described above, a set of 'forward detectors' are also 
used in the present analysis. The copper / silicon Plug calorimeter, the Forward Muon Detector 
(FMD) and the Proton Remnant Tagger (PRT) are sensitive to hadronic activity at large pseudo- 
rapidity, near to the outgoing proton beam. The Plug enables energy measurements to be made 
in the pseudorapidity range 3.5 < i] < 5.5. The FMD is a series of drift chambers covering 
the range 1.9 < r\ < 3.7. Primary particles produced at larger r\ are often detected indirectly 
in the FMD if they undergo a secondary scattering with the beam-pipe or other material. The 
PRT, a set of scintillators surrounding the beam pipe at z = 26 m, detects charged particles in 
the region 6.0 < rj < 7.5. 

The ep luminosity is determined by measuring the rate of the QED Bremsstrahlung process 
ep — > ep7 using a photon calorimeter close to the backward beam pipe at z = —103 m. A much 
more detailed description of the HI apparatus can be found in [43,44]. 
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3.2 Data Samples 



Different event samples are used for different Q 2 ranges of the measurement, as summarised 
in table 1. For the interval 3 < Q 2 < 13.5 GeV 2 , a 'minimum bias' sample ('1997 MB') 
corresponding to an integrated luminosity of 2.0 pb _1 is used, which was recorded during a 
dedicated data taking period in 1997 with unbiased triggers. For intermediate photon virtualities 
(13.5 < Q 2 < 105 GeV 2 ), data taken throughout 1997 are used ('1997 all'), corresponding to 
an integrated luminosity of 10.6 pb -1 . The kinematic range Q 2 > 133 GeV 2 is covered by a 
sample ('1999-2000') corresponding to 61.6 pb~\ taken in the years 1999 and 2000. 

For all three samples, HERA collided positrons with protons, the positron beam energy 
being E e = 27.5 GeV in each case. The proton beam energy was E p = 820 GeV in 1997 and 
E p = 920 GeV in 1999 and 2000, leading to ep centre of mass energies of y/s = 301 GeV and 
y/s = 319 GeV, respectively. The '1997 MB' and '1997 all' samples are used to study neutral 
current interactions with the scattered electron 2 detected in the SpaCal calorimeter. More details 
of the analysis of diffractive data with SpaCal electrons may be found in [45]. The '1999-2000' 
sample is used for the study of both neutral and charged current interactions. In the neutral 
current case, the scattered electron is detected in the LAr calorimeter. These measurements are 
further described in [46]. 



Q 2 range 


Data Set 


Proton Energy E p 


Luminosity 


3 < Q 2 < 13.5 GeV 2 
13.5 < Q 2 < 105 GeV 2 
Q 2 > 133 GeV 2 


1997 MB 
1997 all 
1999-2000 


820 GeV 
820 GeV 
920 GeV 


2.0 pb" 1 
10.6 pb" 1 
61.6 pb" 1 



Table 1: Summary of the data samples used in the analysis. 



3.3 Selection and Reconstruction of DIS Events 

The trigger conditions, detector alignment and calibration, and inclusive DIS selection criteria 
are very similar to those used in the analogous fully inclusive HI analyses at low [47] and 
high [41] Q 2 . The selection criteria are summarised below. 

Neutral current DIS events are triggered by the detection of an energetic electromagnetic 
calorimeter cluster attributed to the scattered electron. Inefficient regions of the calorimeters, 
for example due to cracks between modules or poorly performing trigger cells, are not included 
in the analysis. The trigger efficiency is then close to 100% for the ranges in electron energy 
considered here, namely E' e > 6.5 GeV for electrons detected in the SpaCal and E' e > 11 GeV 
for LAr electrons. To suppress photoproduction background, in which the scattered electron 
escapes undetected in the backward direction and a hadron fakes the electron signature, cuts are 
applied on the lateral extent and isolation of the cluster forming the electron candidate and its 
containment within the electromagnetic part of the calorimeter. 

2 The scattered positron is referred to as an electron throughout this paper. 
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An event vertex, reconstructed by the central or forward tracker, is required within 30 (35) cm 
of the nominal interaction point for the SpaCal (LAr) electron samples. To suppress background 
where a photon fakes the scattered electron, a charged particle track segment must be associ- 
ated to the electron candidate. This track is reconstructed in the BDC for SpaCal electron 
candidates and in the central tracking system for LAr electron candidates. The polar angle 9' e 
of the scattered electron is calculated from the interaction vertex and the BDC track (the LAr 
cluster) for SpaCal (LAr) electrons. In order to ensure good acceptance for the electron in 
the calorimeters and associated trackers, the scattered electron polar angle range considered is 
156° < 9' e < 176.5° for the SpaCal electron samples and 9' e < 153° for the LAr electron sample. 

Hadronic final state objects are reconstructed from the LAr and SpaCal calorimeters and 
the central tracking system using an energy flow algorithm which combines charged particle 
tracks with calorimeter deposits without double counting [48]. Isolated low energy calorimeter 
deposits are classified as noise and are rejected from the analysis. To further suppress back- 
grounds, a minimum of two remaining hadronic final state objects is demanded. Consistency is 
required between the variable y, as reconstructed according to 

y e = 1 - E' e /E e ■ sin 2 (^/2) (Electron method) , 

y h = (E—p z ) h /2E e (Hadron method) , 

y d = tan(7/2)/ [tan(^/2) + tan(7/2)] (Double angle method) , 

using, respectively, the scattered electron only, the hadronic final state 3 only [49], and the elec- 
tron and hadronic final state polar angles 9' e and 7 [50]. The criteria \y e — y h \ < 0.25 and 
\Ue — Ud\ < 0.25 remove badly reconstructed events, further suppress photoproduction back- 
ground and reduce QED radiative corrections due to photon emission from the initial state 
positron. 

The variables y, Q 2 and x are reconstructed by combining information from the scattered 
electron and the hadronic final state using the method introduced in [5]: 

y = y e + y d (i-y d ); Q = tan2(g , /2) ; * = -■ (6) 

This reconstruction method interpolates between the electron method at large y where it has the 
best performance and the double angle method at low y. For diffractive events, where the X 
component of the hadronic final state is well contained in the central detectors, the polar angle 
7 is well measured and the double angle method has an improved resolution compared with that 
for non-diffractive events. This method yields a resolution of 5 — 15% in y, improving with 
increasing y. The resolution in Q 2 is around 3%. 

The principal selection criterion for the charged current sample is a large missing transverse 
momentum, corresponding to the unobserved final state neutrino. This is identified at the trigger 
level mainly using the Liquid Argon calorimeter. The trigger efficiency exceeds 60% throughout 
the kinematic range studied here. For the final selection, the missing transverse momentum must 
exceed 12 GeV. A reconstructed event vertex is required as for the neutral current case. To 

3 Here and in the following, the four- vector of the reconstructed hadronic final state is denoted (E,p x ,p y ,p z )f l . 
In diffractive events, ep — > eXY, where the leading baryonic system Y is not observed, this corresponds to the 
four- vector of the X system. 
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suppress backgrounds, further selection criteria are applied on the event topology, as described 
in [41]. The inclusive kinematic variables are reconstructed using the final state hadrons [49]. 

The selection y > 0.04 is applied to all data samples to ensure reasonable containment of 
the hadronic final state in the central detectors. For the neutral current LAr electron data, the 
sample is restricted to y e < 0.63 for Q 2 < 890 GeV 2 and to y e < 0.9 for Q 2 > 890 GeV 2 , 
which suppresses photoproduction background. 



3.4 Selection and Reconstruction of Diffractive Events 

Diffractive events are selected on the basis of a large rapidity gap, separating the leading bary- 
onic system Y from the system X. The rapidity gap is inferred from the absence of activity in 
detectors sensitive to forward energy flow. The pseudorapidity of the most forward energy de- 
posit above a noise threshold of 400 MeV in the LAr calorimeter must satisfy rj max < 3.3. This 
requirement ensures that the forward extent of the X system lies within the acceptance range 
of the main detector components and thus that M x can be reconstructed reliably. There must 
also be no activity above noise thresholds in the Plug, FMD and PRT detectors. Studies of the 
correlations between the activity levels in the different forward detector components show that 
the LAr, Plug, FMD and PRT requirements have rejection efficiencies for events with no large 
rapidity gap of around 95%, 80%, 80% and 30%, respectively. These efficiencies are well de- 
scribed by the simulations of inclusive DIS used in the analysis (section 3.6) and the combined 
efficiency for the rejection of events with hadronic activity in the range 3.3 < rj < 7.5 is close 
to 100%. Corrections of around 5%, evaluated using randomly triggered events, are applied to 
account for the component of the diffractive signal which is rejected due to electronic noise, 
synchrotron radiation or other effects which fake activity in the forward detectors. 

The diffractive event kinematics are reconstructed using the mass of the system X, obtained 
from 



M x = ^(E 2 -p 2 x -p 2 y -p 2 z ) h 



- ■ (7) 

Vh 



Neglecting the transverse momentum of the hadrons, this method of reconstructing M x reduces 



at large y to a measurement of y 2(E + p z )h (E e — E' e sin 2 ^f) , thus improving the resolution 
where losses of hadrons in the backward (—z) direction become significant. The resolution in 
M x varies between 20% and 30% in the measured kinematic range. The diffractive variables (3 
and x p are obtained from 



« Q 2 x m 

fi = vTJq' x ' = ~r (8) 

The sensitivity of the measurement to variations in the details of the selection and recon- 
struction has been tested in detail. For example, there is no significant change in the results 
when one or more of the forward detectors is not used in the measurement. The final charged 
current sample is scanned visually and all events are consistent with production via ep scattering 
in the interaction region. 
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3.5 Measurement Ranges in M Y and t 



The large rapidity gap selection yields a sample which is dominated by the single dissociation 
process ep — > eXp, with the proton transverse momentum p t (p), and hence \t\ ~ p 2 (p), rela- 
tively small. However, there is an admixture of proton dissociative events, ep — > eW, where 
the proton dissociation system has a small mass M Y . The ranges of sensitivity of the measure- 
ment in M Y and t are determined by the acceptances of the forward detectors which are used to 
identify the large rapidity gap (section 3.4). 

In order to keep the uncertainties arising from proton dissociation small and to ease com- 
parisons with previous data [5, 33, 34], the measurement is integrated over the region 

M Y < 1.6 GeV , \t\ < 1 GeV 2 . (9) 

The correction factors applied to account for the net migrations about these limits are deter- 
mined by evaluation of the forward detector response to elastic proton and proton dissociative 
processes 4 , using the DIFFVM [51] simulation. The ratio of the full generated proton dissocia- 
tion cross section to the generated elastic cross section is taken to be 1:1. Proton dissociation is 
simulated using an approximate da/dM 2 oc 1/M 2 dependence, with explicit simulations of the 
most important resonances at low M Y [52]. The t dependence for proton dissociation follows 
da/dt oc e Bpot with a slope parameter B PD = 1 GeV -2 . The uncertainties are evaluated by 
varying the details of this simulation as described in section 3.7. The resulting correction factors 
are -8.2 ± 5.8% ('1997 MB'), -8.6 ± 5.8% ('1997 all') and -12.0 ± 7.4% ('1999-2000'). 

Comparison of the current data with a similar measurement in which the leading proton is 
directly measured [6] yields a ratio of cross sections for M Y < 1.6 GeV and M Y = m p of 
1.23±0.03 (stat.) ±0.16 (syst.) which is consistent with the DIFFVM prediction of 1.15l{j;£f. 
Neither the comparison in [6], nor further studies sensitive to the larger M Y region [45], show 
any evidence for a dependence of the ratio of proton dissociation to elastic cross sections on any 
of the measured kinematic variables in the region under study. 



3.6 Simulations and Corrections to the Data 

Corrections for detector inefficiencies, acceptances and migrations between x F , (3 and Q 2 mea- 
surement intervals are evaluated using a Monte Carlo simulation which combines several mod- 
els. The RAPGAP [53] event generator simulates the processes e + p — > e + Xp and e + p — > v e Xp 
with x F < 0.15, assuming proton vertex factorisation. Both pomeron and sub-leading ex- 
changes are included, based on the DPDFs from a leading order QCD fit to previous HI data 
('fit 2' in [5]). The parton densities are evolved using Q 2 as a scale and are convoluted with 
leading order QCD matrix elements. Higher order QCD radiation is modelled using initial and 
final state parton showers in the leading \og(Q 2 ) approximation [54]. Hadronisation is simu- 
lated using the Lund string model [55] as implemented in JETSET [56]. QED radiative effects, 
including virtual loop corrections, are taken into account via an interface to the HERACLES 
program [57]. Small weighting factors are applied to the neutral current simulation to ensure 

4 Only proton dissociation to low M Y states is considered in this procedure. Proton dissociation with M Y > 
5 GeV is simulated using the inclusive Monte Carlo models as described in section 3.6. 
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that the t dependence matches that measured in [6] and to optimise the description of the current 
data. The DIFFVM model [51] is used to simulate the exclusive production of the p, u, and 
J/ip vector mesons, which contribute significantly in the SpaCal electron samples at small M x 
(high/5). 

Due to the small inefficiency in the rejection of events with forward hadronic activity using 
the forward detectors, small contributions in the selected data arise from the regions x p > 0.15 
and M Y > 5 GeV. These backgrounds are subtracted using a simulation based on the DJANGO 
[58] Monte Carlo model of inclusive DIS for the SpaCal electron sample and the non-diffractive 
RAPGAP simulation [53] for the LAr electron and charged current samples. The cross sections 
in the inclusive simulations are obtained from QCD fits to recent HI DIS data [47, 59]. 

Residual photoproduction background, which is sizeable only at the highest y values, is 
subtracted on a statistical basis. Its contribution to the data is evaluated using the PHOJET [60] 
model for the SpaCal electron data and the PYTHIA [56] Monte Carlo model for the charged 
current sample. For the neutral current LAr data, the background is evaluated from the fraction 
of reconstructed events for which the detected lepton candidate has the opposite charge to the 
beam lepton, under the assumption that the background is charge symmetric [41]. The small 
backgrounds near (3=1 from QED-Compton scattering (ep — > ep'j) and from di-lepton pro- 
duction via photon-photon fusion (ep — > epe + e~) are subtracted using the COMPTON [61] 
and LPAIR [62] Monte Carlo simulations, respectively. In the charged current measurement, 
a further small background from the production of real electroweak gauge bosons is simulated 
using the EPVEC [63] Monte Carlo model. The normalisations and kinematic dependences of 
each of the background simulations have been checked using dedicated alternative selections 
designed to enhance the corresponding background. 

3.7 Systematic Uncertainties 

A detailed systematic error analysis has been performed, in which the sensitivity of the mea- 
surements to variations in the efficiencies and energy scales of the detector components and 
to the details of the correction procedure is tested. For the neutral current measurement, the 
systematic error sources leading to uncertainties which are correlated between data points are 
listed below. 

• The uncertainty on the SpaCal electromagnetic energy scale varies from 2.4% at E e = 
6.5 GeV to 0.5% at E e = 27.5 GeV [47]. For electrons detected in the LAr calorimeter, 
the energy scale is known to a precision varying between 1% and 1.5%, depending on 
9' e [41]. The uncertainties in the relative alignment of the different detector components 
are reflected in possible biases in the electron polar angle measurement at the level of 
0.5 mrad for the SpaCal data [47] and between 1 mrad and 2 mrad, depending on 9' e , for 
the LAr data [41]. 

• The hadronic energy scale of the LAr calorimeter is known to 2% for all samples stud- 
ied. The uncertainty on the hadronic energy scale of the SpaCal is 7% and that on the 
contribution to the hadronic energy measurement from charged particle tracks is 3%. 
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• Imperfect treatment of calorimeter noise can result in a bias in the reconstruction of M x . 
The corresponding uncertainty is evaluated by varying the amount of calorimeter en- 
ergy classified as noise by 10%. This level of precision is determined by comparing the 
calorimeter noise subtracted in the data with the Monte Carlo model, which includes a 
simulation of noise based on randomly triggered events. 

• The efficiency with which the FMD registers activity when there is hadronic energy flow 
in its acceptance region is varied by 5% in the simulation. For the PRT, this efficiency is 
varied by 20%. The Plug energy scale is varied by 30%. These levels of uncertainty are 
obtained by comparison of the present data with the Monte Carlo simulation for samples 
in which forward detector activity is required to be present rather than absent. 

• The model dependence of the acceptance and migration corrections and background sub- 
tractions is estimated by varying the details of the Monte Carlo simulation within the 
limits permitted by the present data. In the RAPGAP simulation of diffraction for the 
SpaCal electron data, the x p distribution is reweighted by (1 / x P ) ±om , the (3 distribution 
by (3 ±om and (l-(3) ±om , the t distribution by e ±l and the Q 2 distribution by (logQ 2 ) ±02 . 
The same systematic shifts are applied for the LAr electron data, except that the powers 
by which the x p , f3 and 1 — f3 distributions are reweighted are increased from 0.05 to 
0.1, reflecting the weaker constraints on those distributions from the data at high Q 2 . 
The normalisation of the sub-leading meson exchange in RAPGAP is varied by ±25% 
and that of the vector meson production simulation (DIFFVM) is varied by ±50%. The 
uncertainty in the background from high x p or M Y , as simulated by the DJANGO and 
inclusive RAPGAP Monte Carlo models, is taken to be 100%. Appropriate variations are 
also made in the normalisations of the photoproduction, QED-Compton scattering and 
lepton pair production background simulations. 

Several further uncertainties, listed below, affect all data points in an identical manner and 
are thus considered as normalisation uncertainties. 

• The uncertainty on the factor accounting for smearing about the M Y and t boundaries 
of the measurement (section 3.5) is 5.8% for the data taken in 1997 and 7.4% for the 
1999-2000 data. The dominant contribution to this uncertainty arises from variations in 
the assumed ratio of proton dissociation to elastic proton cross sections in the range 0.5 
to 2.0, which is determined from studies with alternative forward detector requirements. 
Smaller contributions arise from reweighting the M Y and t distributions of the proton 
dissociation simulation and propagating the uncertainties on the FMD, PRT and Plug 
detectors. 

• The uncertainty arising from the luminosity measurement is 1.5% for all samples. 

• The correction factor for diffractive events rejected due to noise fluctuations in the for- 
ward detectors is varied by 25%, which corresponds to the r.m.s. variation over the dif- 
ferent fills of HERA. This leads to normalisation uncertainties at the 1% level, varying 
slightly between the different data sets. 
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A final class of systematic errors leads to uncertainties which are not taken to be correlated 
between data points. 

• The calculated acceptance of the r) max cut depends on the modelling of the hadronic fi- 
nal state topology. The associated uncertainty is estimated from the effect of using an 
alternative model for higher order QCD processes (the colour dipole approach [64] as 
implemented in ARIADNE [65] in place of parton showers). This results in an uncer- 
tainty which depends to good approximation on x p only and varies between 1.2% at 
x p = 0.0003 and 11% at x p = 0.03. 

• The uncertainty in the trigger efficiency is 1.0% in the SpaCal data [47] and 0.3% for LAr 
data [41]. 

• Uncertainties of 0.5% arise in all neutral current data sets due to the uncertainty in the 
efficiency of the track-link requirements for the electron candidate [41,47]. 

The total systematic uncertainty on each data point is formed by adding the individual con- 
tributions in quadrature. Away from the boundaries of the kinematic region studied in the neu- 
tral current measurement, the systematic error varying from point to point is around 5%, with 
no single source of uncertainty dominating. The systematic error increases to typically 15% 
at the largest x p = 0.03, where the contribution from the modelling of the acceptance of the 
rj max requirement becomes important. At small M x values, the calorimeter noise uncertainty 
becomes the largest. These point-to-point systematic uncertainties are to be compared with nor- 
malisation uncertainties of 6.2% and 7.6% for the SpaCal and LAr electron data, respectively, 
and statistical errors of between 5% and 20%. 

For the charged current measurement, the statistical uncertainties are dominant. The sys- 
tematic error treatment is similar to that used in the neutral current case, except that the model 
variations are increased in light of the lack of previous data to constrain the kinematic depen- 
dences. The largest systematic uncertainties arise from these model uncertainties, in particular 
from that on the sub-leading meson exchange contribution, and from the modelling of the ac- 
ceptance of the r] mSLX requirement. 



4 The Diffractive Neutral Current Cross Section 

In order to obtain the reduced neutral current cross section defined in equation 4, the data are 
corrected to fixed values of Q 2 , (3 and x p , the influence of the finite bin sizes being evaluated 
using a parameterisation of the QCD fit described in section 5. The measurements are quoted 
at the Born level after the corrections for QED radiative effects described in section 3.6. For all 
data points shown, the acceptance, bin purity and bin stability 5 exceed 30 %. The results are 
given in numerical form in tables 4 (SpaCal electron data) and 5 (LAr electron data) and are 
shown graphically as described in sections 4.1 and 4.2. They can also be found at [66]. 

5 Purity (stability) is evaluated using the Monte Carlo simulation and is denned as the fraction of events recon- 
structed (generated) in a measurement bin which are also generated (reconstructed) in that bin. 
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4.1 Dependences on (3 and Q 2 



The QCD properties and structure of the diffractive interaction are most easily interpreted from 
the dependences on x (or equivalently f3) and Q 2 , with x p fixed. According to [20], DPDFs 
can then be defined for each fixed x p value, independently of the validity of proton vertex 
factorisation. A binning scheme is therefore adopted with fixed x p , x and Q 2 values. In order 
to minimise the statistical uncertainties and the systematic uncertainties associated with the 
reconstruction of M x , relatively large x p intervals are chosen. The binning in the much better 
resolved variables, x and Q 2 , is chosen to match previous inclusive measurements [47]. 

The (3 and Q 2 dependences of the reduced cross section, multiplied by x p , are shown in 
figures 2-6 at fixed values of x p = 0.0003, 0.001, 0.003, 0.01 and 0.03, respectively. The 
results from different x p values complement one another in their j3 and Q 2 ranges, though there 
is also considerable overlap between the coverage at different x p values. For each x p value 
considered, the data exhibit similar (3 and Q 2 dependences. 

As can be seen in figures 2a-6a, the cross section remains large up to the highest accessed 
values of (3 (i.e. where x — > x F ) at fixed x p and Q 2 . This behaviour is in marked contrast 
to that of hadron structure functions, but is qualitatively similar to that of the photon [67]. 
The Q 2 dependence of the data is shown for fixed x p and (3 in figures 2b-6b. The reduced 
cross section increases with Q 2 throughout most of the kinematic range, up to large f3 ~ 0.5. 
These positive scaling violations confirm earlier observations [5] and contrast with the case of 
inclusive scattering from hadrons, for which the cross section at fixed x falls with increasing Q 2 
forx> 0.1 [47]. 

The data in figures 2-6 are compared 6 with the results of the 'HI 2006 DPDF Fit A de- 
scribed in section 5. This fit assumes proton vertex factorisation and includes a sub-leading 
exchange contribution, which is important at low (3 and large x p as shown in figures 5a and 6a. 
It is clear from the good overall description that the data are broadly consistent with such a 
model. The (3 dependence of then directly reflects the quark structure of the diffractive 
exchange with each quark flavour weighted by its squared electric charge, whilst the measured 
Q 2 dependence is sensitive to the diffractive gluon density. 

The Q 2 dependence is quantified by fitting the data at fixed x p and (3 to the form 

v^\x p ,Q 2 ,[3) = a D ([3,x p ) + b D ([3,x p )\nQ 2 , (10) 

such that bn{(3,x p ) = da®^ /d\nQ 2 is the first logarithmic Q 2 derivative of the reduced 

cross section. This observable has been used previously to discriminate between different mod- 
els of diffractive DIS [68]. Equation 10 is fitted to data with 0.001 < x p < 0.03 from each 
(3 value if there are a minimum of three available data points 7 . The resulting \nQ 2 derivatives 
are shown in figure 7a, after dividing b D ((3, x p ) by the factor fjp/ p (x p ), defined in equation 14, 
which is used to parameterise the x p dependence so that the results from different x p values 



6 The curves shown correspond to E p = 820 GeV. The predictions for E p = 920 GcV differ slightly at the 
lowest (3 values, due to the influence of F^ 3 ' '. 

7 If only data with Q 2 > 8.5 GeV 2 are included, as in the fit described in section 5, the changes to the logarith- 
mic Q 2 derivatives are small and the conclusions are unaffected. 
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can be compared in normalisation as well as in shape. Although the logarithmic derivatives at 
different x p values probe different Q 2 regions, they are remarkably similar when viewed as a 
function of (3. This confirms the applicability of the proton vertex factorisation framework to 
the description of the current data. The lack of any significant change in behaviour at large 
x p indicates that the derivatives are not significantly affected by the presence of sub-leading 
exchange contributions. 

According to the DGLAP evolution equations, the InQ 2 derivative of F 2 D< ^ contains con- 
tributions due to the splittings g — > qq and q — > qg, convoluted with the diffractive gluon and 
quark densities, respectively. In figure 7b, the derivatives at x p = 0.01 are shown together 
with the decomposition into these two contributions according to the 'HI 2006 DPDF Fit A' 
described in section 5. The curves correspond to Q 2 values which vary in order to match the 
average In Q 2 of the data which are fitted at each (3 value. The theoretical calculation is in good 
agreement with the data. The In Q 2 derivative is determined almost entirely by the diffractive 
gluon density up to (3 ~ 0.3. The large positive InQ 2 derivatives in this region can thus be 
attributed to a large gluonic component in the DPDFs. For (3 > 0.3, the contribution to the Q 2 
evolution from quark splittings becomes increasingly important and the derivatives become less 
sensitive to the gluon density. The fall in the derivative as (3 — > 1 then arises dominantly from 
gluon radiation, q — > qg, shifting quarks from higher to lower f3 with increasing Q 2 . 

Considering diffractive DIS in terms of the elastic scattering from the proton of colour 
dipoles formed by partonic fluctuations of the exchanged virtual photon [69] provides a comple- 
mentary framework in which to describe diffractive DIS. Modelling the exchange by two gluons 
in a net colour singlet configuration [70], the data at low and moderate (3 are then described in 
terms of qq and qqg fluctuations of transversely polarised photons, whilst the high (3 region 
contains a Q 2 -suppressed non-leading twist contribution from qq fluctuations of longitudinally 
polarised photons [17, 18,71]. A significant contribution from this term could explain the ten- 
dency of the data to rise with increasing f3 for f3 > 0.1 at low Q 2 (figure 2a- 4a). However, the 
high (3 data in figures 2b-6b are adequately described by a single logarithmic dependence on Q 2 
and do not require a sum of leading and Q 2 -suppressed terms. 

4.2 Dependence on x p and Comparisons with other Data 

For comparison with previous measurements of diffractive DIS and for a more detailed study of 
the x p dependence at fixed (3 and Q 2 , the reduced cross section is also measured using a binning 
scheme with fixed (3, Q 2 and x = (3 ■ x p , as shown in figure 8. The data are multiplied by x p 
for visibility and are compared with the results of the QCD fit described in section 5, which is 
also in good agreement with the data obtained using this binning scheme. The x p dependence 
is roughly flat for all (3 and Q 2 values, implying that the reduced cross section approximately 
follows a <7r^ oc l/x p dependence. However, when viewed in detail, there are clear deviations 
from this behaviour. The variations in the x p dependence as (3 changes are as expected from 
the interplay between the leading pomeron and a sub-leading trajectory exchange, as discussed 
in section 5.4. 

In figure 9, the measured reduced cross section is compared with results obtained by di- 
rect measurement of the final state proton using the HI Forward Proton Spectrometer ('HI 
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(FPS)') [6] and the ZEUS Leading Proton Spectrometer ('ZEUS (LPS)') [10]. In figure 10, a 
comparison is made with ZEUS data obtained by decomposition of the inclusive In M 2 distri- 
bution into diffractive and non-diffractive components ('ZEUS (M x )') [11]. Shifts, evaluated 
using the DPDF fit described in section 5, are applied to these data in order to transport them 
to the p and Q 2 values of the present measurement. Since no uncertainties are ascribed to this 
procedure, only data points for which the shifts are small and relatively insensitive to the choice 
of DPDFs are shown. The different contributions from proton dissociation in the different data 
sets are accounted for by the application of global factors of 1.23 to the 'HI (FPS)' and 'ZEUS 
(LPS)' data and of 0.86 to the 'ZEUS (M x )' data. The former factor corresponds to the mea- 
sured ratio of cross sections for M Y < 1.6 GeV and M Y = m p , for which the uncertainty is 
0.16 (see [6] and section 3.5). The latter factor contains an additional contribution of 0.7, cor- 
responding to the ratio of cross sections for M Y = m p and M Y < 2.3 GeV according to [1 1]. 

There is broad agreement between all of the data sets on the general behaviour of the diffrac- 
tive cross section. The 'HI (FPS)' and 'ZEUS (LPS)' data are compatible with the present mea- 
surement throughout the kinematic range available for comparison. A more detailed comparison 
with the 'HI (FPS)' data is presented in [6]. The 'ZEUS (M x )' data are in good agreement with 
the present measurement in some regions, for example at low (3 and low Q 2 . However, there are 
disagreements in the low Q 2 , high (3 region and in the high Q 2 , low (3 region, which correspond 
to low and high values of M x = >/Q 2 (1/(3 — 1), respectively. The high M x discrepancy gives 
rise to a weaker Q 2 dependence of the 'ZEUS (M x )' cross section at low f3 than is the case for 
HI [72]. 



5 QCD Analysis and Diffractive Parton Distributions 

The high precision and large kinematic range of the diffractive cross section data presented in 
this paper allow detailed tests of the factorisation properties of diffractive DIS and the extraction 
of DPDFs, which may be used to predict cross sections for other diffractive processes at HERA 
and elsewhere. 

5.1 Theoretical Framework 

QCD hard scattering collinear factorisation, when applied to diffractive DIS [20], implies that 
the cross section for the process ep — > eXY can be written in terms of convolutions of partonic 
cross sections <j ei (x, Q 2 ) with DPDFs ff as 

da^ eXY (x,Q 2 ,x F ,t) =J2 tf(x,Q 2 ,x F ,t) ® da™(x,Q 2 ) . (11) 

i 

The partonic cross sections are the same as those for inclusive DIS. The DPDFs represent prob- 
ability distributions for partons i in the proton under the constraint that the proton is scattered 
to a particular system Y with a specified four-momentum. They are not known from first prin- 
ciples, but can be determined from fits to the data using the DGLAP [39] evolution equations. 
The factorisation formula in equation 11 is valid for sufficiently large Q 2 and fixed x p , t and 
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system Y. It also applies to any cross section which is integrated over a fixed range in M Y and 
t and may thus be applied to the present data with M Y < 1.6 GeV and \t\ < 1 GeV 2 . 

Due to kinematic constraints, it is not possible to access the full range of x and Q 2 using data 
from only one value of x p . A parameterisation of the x p dependence of the DPDFs, for which 
there is no clear procedure in QCD, is therefore necessary. The proton vertex factorisation 
framework is adopted here, such that the DPDFs are factorised into a term depending only on 
x p and t and a term depending only on x (or (3) and Q 2 : 

f?(x,Q 2 ,x p ,t) = fjp /p (x p ,t) ■ f0 = x/x p ,Q 2 ) . (12) 

This is equivalent to treating the diffractive exchange as a 'pomeron' with a partonic struc- 
ture given by the parton distributions fi((3,Q 2 ), the variable f3 corresponding to the fraction 
of the pomeron longitudinal momentum carried by the struck quark. The 'pomeron flux factor' 
f]p/ p (x p , t) represents the probability that a pomeron with particular values of x p and t couples 
to the proton. In the low x p region where sub-leading exchange contributions are negligible, the 
data presented in section 4 are consistent with factorisation of the x p dependence, as required 
for the proton vertex factorisation expressed in equation 12. There is also no evidence from the 
investigations in [6, 10] for any change in the t or M Y dependences as either j3 or Q 2 vary in the 
range relevant to the present analysis. 



5.2 Fit Procedure 

To determine the DPDFs, fits are made to the reduced cross section data as presented in fig- 
ures 2-6, for which (3 < 0.8. In order to avoid regions which are most likely to be influenced 
by higher twist contributions or other problems with the chosen theoretical framework, only 
data with M x > 2 GeV are included in the fit and the region Q 2 < 8.5 GeV 2 is excluded, as 
explained in section 5.3. The total number of fitted data points is 190. 

Input parameters describing the DPDFs at a starting scale Ql for QCD evolution are ad- 
justed to obtain the best description of the data after NLO DGLAP [40] evolution 8 to Q 2 > Ql 
and convolution of the DPDFs with coefficient functions. The fit is performed in the MS renor- 
malisation scheme [74] with charm and beauty quarks treated as massive, appearing via boson 
gluon fusion-type processes up to order a 2 [75]. The heavy quark masses are set to world aver- 
age values (see table 2). The strong coupling is set via Aq^, d = 399 ± 37 MeV for 3 flavours, 
which corresponds [76] to the world average af\M^) for five flavours (table 2). The effects 
of are considered through its relation to the quark and gluon densities at NLO. Since the 
relative normalisations of the diffractive parton densities and the flux factor in equation 12 are 
arbitrary, no momentum sum rule is imposed. 

The DPDFs are modelled in terms of a light flavour singlet distribution consisting of 
u, d and s quarks and anti-quarks with u = d = s = u = d = s, and a gluon distribution 
g(z). Here, z is the longitudinal momentum fraction of the parton entering the hard sub-process 
with respect to the diffractive exchange, such that z = (3 for the lowest order quark-parton 

8 Direct pomeron to parton splitting functions, leading to an inhomogeneous term in the DGLAP evolution 
equations [32, 73], are not considered. The presence of such a term would lead to a reduced gluon density. 
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model process, whereas < (3 < z for higher order processes. The quark singlet and gluon 
distributions are parameterised at Q\ using a similar approach to that commonly applied to 
hadronic parton densities [77-79], such that the most general form is 

zf i (z,Q 2 )=A i z Bi (l-z) Ci . (13) 

The exact choices of terms included for the quark singlet and gluon densities and the param- 
eterisation scale Ql are determined through a systematic investigation of the parameter space, 
as described in section 5.3. The DPDFs as defined in equation 13 are multiplied by a term 

0.01 

e in order to ensure that they vanish at z = 1, as required for the evolution equations to be 
solvable. The parameters C q and C g thus have the freedom to take negative as well as positive 
values. Modifying the argument of the exponential term within reasonable limits has no visible 
influence on the fit quality or the extracted DPDFs in the range of the measurement. 

The x p dependence is parameterised using a flux factor motivated by Regge theory, 

fjp/p(x F ,t) = Ajp ■ 2Qf(f) _ 1 , (14) 

where the pomeron trajectory is assumed to be linear, ajp(t) = q;jp(0) + a'jpt, and the pa- 
rameters Bjp and a' p and their uncertainties are obtained from fits to HI FPS data [6]. The 
values of these and other parameters which are fixed in the fits are summarised in table 2. Fol- 
lowing the convention of [5], the value of the normalisation parameter Ajp is chosen such that 
x p ' ft™™ fjP/p dt = 1 at = 0.003, where |i m i n | — m p x F / (1 — x F ) is the minimum kine- 
matically accessible value of \t\, m p is the proton mass and \t cut \ = 1.0 GeV 2 is the limit of the 
measurement. 



Parameter 


Value 


Source 


a'jp 
Bp 
«jr(0) 

BjR 

m c 
m b 
af\Ml) 


0.06 ±g;J§ GeV" 2 
5.5 - 2 ° 7 GeV" 2 

0.50 ±0.10 
0.3 t 6 3 GeV- 2 
1.6 ~li GeV- 2 

1.4 ±0.2 GeV 

4.5 ±0.5 GeV 
0.118±0.002 


[6] 

[6] 

[5] 

[6] 

[6] 

[52] 

[52] 

[52] 



Table 2: The values of the fixed parameters and their uncertainties, as used in the QCD fits. 
Since they are strongly anti-correlated when extracted from fits to the FPS data, a'jp and Bjp 
are varied simultaneously to obtain the theoretical errors on the fit results, as are a'jp and Bjp. 
The remaining parameters are varied independently. The theoretical uncertainties on the free 
parameters of the fit also contain a contribution from variation of the parameterisation scale Ql, 
as described in section 5.3. 

To obtain a good description of the data, an additional sub-leading exchange (M) is included, 
which has a lower trajectory intercept than the pomeron and which contributes significantly only 
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at low (3 and large x F . As in [5,6] this contribution is assumed to factorise in the same way as 
the pomeron term, such that equation 12 is modified to 

f?(x,Q\x p ,t) = fjp /p {x p ,t)- f^^Q 2 ) + nn-fn/^x^t)-^^ 2 ). (15) 

The flux factor f M / p takes the form of equation 14, normalised via a parameter Am in the same 
manner as for the pomeron contribution and with fixed parameters or(0), a' M and Bjr obtained 
from other HI measurements (see table 2). The parton densities ff* of the sub-leading ex- 
change are taken from a parameterisation derived from fits to pion structure function data [80]. 
Choosing a different parameterisation [81] does not affect the fit results significantly. 

The free parameters of the fit are the A, B and C parameters which determine the quark 
singlet and gluon distributions (equation 13), together with 0:^(0), which controls the x p de- 
pendence and n M , which controls the normalisation of the sub-leading exchange contribution. 
In order to constrain these parameters, a x 2 function as defined in [47] is minimised. This func- 
tion involves the combined statistical and uncorrelated systematic errors for each data point and 
also takes account of correlations between data points caused by systematic uncertainties by 
allowing variations in each systematic error source at the expense of increases in the x 2 vi- 
able. Ten correlated systematic error parameters are considered for each of the SpaCal and LAr 
electron data sets, one for each of the error sources listed in tables 4 and 5. In this procedure, the 
sources of correlated uncertainty are taken to be common for the ' 1 997 MB ' and ' 1 997 all' data, 
whereas it is assumed that there are no correlations between the uncertainties on the SpaCal and 
the LAr electron data. A further systematic error parameter controls the relative normalisation 
of the LAr electron data set with respect to the SpaCal electron data, for which the uncertainty 
is 9.8% (section 3.7). The central results for the DPDFs and other parameters obtained from the 
fit are not altered significantly if all systematic uncertainties leading to correlations between the 
data points are ignored. 

The statistical and experimental systematic errors on the data points and their correlations 
are propagated [82] to obtain experimental uncertainties on the DPDFs and other fit parameters, 
which correspond to increases in the x 2 variable by one unit. The theoretical error is obtained 
from variations of the assumed parameters as given in table 2, with an additional contribution 
expressing the sensitivity to the choice of DPDF parameterisation, obtained by varying Ql as 
discussed in section 5.3. Since the pomeron flux factor is constrained simultaneously with the 
parton densities, the possible influence of interference between the pomeron and sub-leading 
exchange contributions cannot be assessed. However, in previous similar fits in which ot p (0) 
was extracted separately from the parton densities [5], a F (0) changed by less than 0.01 between 
the cases of no interference and maximum constructive interference. 

5.3 Choices of Fit Parameterisation and Kinematic Range 

In order to optimise the results of the fit, the sensitivity to variations in the details of the parame- 
terisation is investigated. With the small numbers of parameters used to describe the parton den- 
sities, the x 2 values and the results of the fits are sensitive to the choice of the parameterisation 
scale Ql [47], so that its value must be optimised by x 2 minimisation for each parameterisation 
choice. The QQ-optimised results are then compared in order to make the final parameterisa- 
tion choice. To ensure that the results of the fit are not sensitive to the kinematic range of the 



20 



data included in the fit, the sensitivity of this procedure to variations in the kinematic cuts is 
also tested. All parton density parameterisation changes and kinematic range variations lead to 
extracted values of 0^(0) which are within the experimental uncertainties (see also section 5.4). 

The only significant sensitivity to the boundaries of the chosen kinematic range occurs when 
the minimum Q 2 value of the data included in the fit, Q 2 nin , is varied. Whereas the quark 
distribution remains stable within uncertainties for all (5 min choices, the gluon distribution for 
z < 0.5 increases systematically as Q^ in varies between 3.5 GeV 2 and 8.5 GeV 2 , changing by 
about 40% in total. The x 2 P er degree of freedom also improves steadily as Q 2 l[n varies over this 
range. There is no evidence for any further variation in the gluon density for Qmin > 8.5 GeV 2 . 
The lowest Q 2 data are therefore omitted from the fit and Q min = 8.5 GeV 2 is chosen. The 
Qmm dependence is reflected in figures 2-6 as a tendency for the fit result extrapolated to Q 2 < 
8.5 GeV 2 to lie below the data. The dependence of the gluon density on <5 min may indicate 
inadequacies in the adopted formalism at the lowest Q 2 values. 

The fit results are not sensitive to variations in the minimum or maximum (3 values of the 
data included, although the minimum (3 is correlated with Qmm through the kinematics. There 
is similarly no significant change in the fit results when the minimum M x value of the data 
included is increased or when the highest x p data are omitted. 

For the quark singlet distribution, the data require the inclusion of all three parameters A q , 
B q and C q in equation 13. By comparison, the gluon density is weakly constrained by the data, 
which are found to be insensitive to the B g parameter. The gluon density is thus parameterised at 
Ql using only the A g and C g parameters. With this parameterisation, a value of Ql = 1.75 GeV 2 
yields the minimum x 2 value of 158 for 183 degrees of freedom. This fit is referred to as the 
'HI 2006 DPDF Fit A' in the figures. As a measure of the parameterisation uncertainty, Ql 
is varied between 1.15 GeV 2 and 2.05 GeV 2 , for which the x 2 variable increases by one unit. 
The correlated systematic error parameters have a mean close to zero and the largest shift in a 
correlated error source is 1.0 a. The fit shifts the normalisation of the LAr data relative to the 
SpaCal data by 0.45 a. Both the DPDFs and the x 2 P er degree of freedom of this fit can be 
reproduced closely using the approach based on Chebyshev polynomials in [5]. 

As discussed in section 4.1 (figure 7b), the Q 2 dependence of the data at fixed (3 and x p 
determines the gluon density well at low (3. However, as (3 increases the In Q 2 derivative be- 
comes smaller and the fractional error on the gluon density becomes correspondingly larger. At 
the highest (3 values, where the Q 2 evolution is driven by quarks, the Q 2 dependence of ov 3 ^ 
becomes insensitive to the gluon density. The results for the gluon density at large z are thus 
determined principally by the data at lower z coupled with the parameterisation choice. This 
lack of sensitivity is confirmed by repeating the fit with the parameter C g , which determines 
the high z behaviour, set to zero. Apart from the exponential term, the gluon density is then a 
simple constant at the starting scale for evolution, which is chosen to be Ql = 2.5 GeV 2 by x 2 
minimisation. Even with this very simple parameterisation of the gluon density, the x 2 variable 
increases only slightly to x 2 = 164, with 184 degrees of freedom. This fit is referred to as the 
'HI 2006 DPDF Fit B' in the figures. 

5.4 Diffractive Parton Distributions and Effective Pomeron Intercept 

A good description of the data is obtained throughout the fitted range Q 2 > 8.5 GeV 2 , (3 < 0.8 
and M x > 2 GeV by both HI 2006 DPDF Fit A and Fit B. The results of Fit A are compared 
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with the measured reduced cross section in figures 2-6. The results for the fit parameters are 
given in table 3. They can also be found together with the correlation coefficients between the 
parameters at [66]. 



1— (if - UorotYicitai* 

rn r didmeier 


rll r\ 


rll D 


a F {0) 


1.118±0.008 


-i 111 i /~i t~\t~\>~7 

l.llliO.007 


n M 


(1.7±0.4) x lO" 3 


(1.4±0.4) x 


A i 


1.06±0.32 


0.70±0.11 


B q 


2.30±0.36 


1.50±0.12 


c q 


0.57±0.15 


0.45±0.09 


Ar 


0.15±0.03 


0.37±0.02 


c 9 


-0.95±0.20 


(fixed) 



Table 3: The central values of the parameters extracted in the 'HI 2006 DPDF Fit A' and 'B', 
and the corresponding experimental uncertainties. 

The diffractive quark singlet and gluon distributions from Fit A are shown together with 
their uncertainties on a logarithmic z scale in figure 1 1 . In order to illustrate the high z region in 
more detail, they are also shown on a linear z scale in figure 12, where they are compared with 
the results from Fit B. At low Q 2 , both the quark singlet and the gluon densities remain large 
up to the highest z values accessed. The quark singlet distribution is well constrained, with an 
uncertainty of typically 5—10% and good agreement between the results of Fit A and Fit B. 
The gluon distribution has a larger uncertainty of typically 15% at low to moderate z and low 
Q 2 , dominated by the influence of the Ql variation. For z > 0.5, where the sensitivity to the 
gluon density becomes poor, the level of agreement between Fit A and Fit B worsens. 

As shown in figure 13, the fraction of the exchanged momentum carried by gluons inte- 
grated over the range 0.0043 < z < 0.8, corresponding approximately to that of the measure- 
ment, is around 70% throughout the Q 2 range studied, confirming the conclusion from earlier 
work [5, 27]. The integrated gluon fraction is somewhat smaller for Fit B due to the lower 
gluon density at large z values, though the results from the two fits are consistent within the 
uncertainties. 

Fit A yields an effective pomeron intercept of 

a„(0) = 1.118 ± 0.008 (exp.) t° mo (model) , (16) 

where the first error is the full experimental uncertainty and the second expresses the model 
dependence. This model dependence uncertainty arises dominantly from the variation of a' , 
which is strongly positively correlated with 0^(0), such that 0^(0) increases to around 1.15 
if a' is set to 0.25. The intercept has also been shown to be strongly sensitive to the value 
of FjP [83], though with the NLO treatment adopted here, is determined in the fit and no 
additional uncertainty is included. The influence of on the reduced cross section according 
to the fit is shown for x p = 0.003 in figure 4a. It is similar at other values of x p . 
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The extracted a F (0) is slightly higher than the value a F (0) ~ 1.08 expected for the 'soft' 
pomeron [84] describing long distance hadronic interactions. The result is compatible with that 
obtained from ZEUS data for Q 2 < 20 GeV 2 [11]. However, in [1 1], evidence was reported 
for an increase of 0^(0) for Q 2 > 20 GeV 2 . In some models [17, 18], a (3 dependent a p (0) 
has also been suggested. Any such dependence of 0^(0) on Q 2 or (3 implies a breakdown of 
proton vertex factorisation. In order to test for such effects in the present data, the QCD fit 
is repeated with additional free parameters corresponding to independent values of a F (0) in 
different ranges of Q 2 or (3. As can be seen from the results in figure 14, there is no evidence for 
any variation of a F (0) with either variable within the kinematic range of the fit. This remains 
the case if the data with Q 2 < 8.5 GeV 2 are included. This lack of Q 2 dependence of a w (0) 
contrasts with the Q 2 dependent effective pomeron intercept extracted in a Regge approach to 
inclusive small x proton structure function data, as discussed further in section 7. 

The presence of the sub-leading exchange term is required by the data, the \ 2 increasing 
by approximately 40 units if only the pomeron contribution is included. Its relative size is 
expressed through the normalisation parameter n R = [1.7 ± 0.4 (exp.) IqJ (model)] x 10~ 3 
(Fit A), where the dominant uncertainty arises from the correlation with 0:^(0). The sub-leading 
exchange plays a significant role at high x F and low (3, as shown in figures 5a and 6a. It accounts 
for around 30% (10%) of the cross section at x F = 0.03 (0.01) and is negligible at lower x F . 



6 The Diffractive Charged Current Cross Section 

The diffractive charged current process e + p — > P e XY is sensitive to the diffractive d, u, s and 
c densities at large scales. Assuming factorisation, the measurement of the charged current 
cross section thus tests the assumed flavour decomposition of the quark singlet component of 
the DPDFs, which is completely unconstrained by the neutral current data. The total charged 
current cross section integrated over the range Q 2 > 200 GeV 2 , y < 0.9 and x p < 0.05 at 
y/s = 319 GeV is measured to be 

<7*g = 390 ± 120 (stat.) ± 70 (syst.) fb , (17) 

corresponding to 2.2 ± 0.7 (stat.) ± 0.4 (syst.)% of the total charged current cross section [41] 
for the same Q 2 and y ranges, with x < 0.05. 

The measured charged current cross section is compared with the prediction of the '2006 
DPDF Fit A' to the neutral current diffractive DIS data described in section 5.4. The prediction 
is obtained by implementing the DPDFs extracted from the neutral current data in the RAPGAP 
[53] Monte Carlo generator. The light quark-initiated contributions are calculated at lowest 
order and the c — > s contribution is calculated using the 0(a s ) matrix element. Leading log 
(Q 2 ) parton showers are used to approximate higher order QCD radiation. The resulting cross 
section prediction is 500 fb, which is compatible with the measurement. The experimental 
uncertainties on the DPDFs and the theoretical uncertainties detailed in table 2 lead to negligible 
errors on the predictions by comparison with the statistical error on the measured cross section. 

The charged current cross section measurement is shown differentially in i p , (3 and Q 2 
in figures 15a-c, respectively. The numerical values are given in table 6. In all cases, the 



23 



predictions derived from the DPDFs of section 5.4 are in agreement with the measurements. 
The charged current data are thus consistent with the singlet quark distribution assumed in the 
DPDF fit, where all light quark and antiquark densities are taken to be equal, although the large 
statistical uncertainties preclude strong conclusions. The contribution in the model from the 
sub-leading exchange is shown as a dashed line in the figures. It contributes at the 15% level 
for x p > 0.015 and is negligible at lower x p . 



7 Comparison between Diffractive and Inclusive DIS 

In hadronic scattering, close connections have been drawn between the diffractive and the to- 
tal cross sections, for example via the generalisation of the optical theorem to diffractive dis- 
sociation processes [85]. These connections are carried forward into many models of low x 
DIS [14, 15, 18]. Comparing the Q 2 and x dynamics of the diffractive with the inclusive cross 
section is therefore a powerful means of developing our understanding of high energy QCD, 
comparing the properties of the DPDFs with their inclusive counterparts and testing models. 

The evolution of the diffractive reduced cross section with Q 2 is compared with that of the 
inclusive DIS reduced cross section ay by forming the ratio 



(18) 



a r (x, Q 2 ) 

at fixed x and x p , using parameterisations of the ay data from 9 [41,47]. This ratio, which 
was also studied in [10], is shown multiplied by x p in figures 16 and 17, as a function of Q 2 
for all measured x F and x = (3x p values. In order to compare the Q 2 dependences of the 
diffractive and the inclusive cross sections quantitatively, the logarithmic derivative of their ra- 
tio, b R (x, xp) = dl ®Q 2 (cy*^ / (T r )x, x p , is extracted through fits of a similar form to equation 10, 
whereby 



OY {% ' p , X, Q j 



(x,Q 2 ) 



= a R (x,x p ) + b R (x,x p ) lnQ 2 . (19) 



The fits are overlayed on the data in figures 16 and 17. The resulting values of 6r are shown in 
figure 18, where they are divided by the flux factor fjp/ p (x p ) (equation 14), to allow compar- 
isons between different x p values. Since the dominant uncertainties arise from the diffractive 
data, the statistical fluctuations in figure 18 reproduce those of figure 7a. 

The ratio of the diffractive to the inclusive cross section is remarkably flat as a function of 
Q 2 for most x and x p values, such that the In Q 2 derivative of the ratio is consistent with zero. 
At the highest (3 > 0.3, where x approaches x p , the ratio falls with increasing Q 2 and the In Q 2 
derivative becomes negative. This occurs in a manner which depends to good approximation on 
P only; at fixed j3, there is no significant dependence of the logarithmic derivative on x p . 

9 The inclusive reduced cross section is denoted <tnc in [41]. 
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The compatibility of the In Q 2 derivative of the ratio with zero over much of the kinematic 
range implies that l/oy 3 ^ • da®^ /d\nQ 2 ~ 1/ay • d<j r /dhiQ 2 . Whereas the diffractive and 
inclusive reduced cross sections are closely related to their respective quark densities, the In Q 2 
derivatives are approximately proportional to the relevant gluon densities in regions where the 
Q 2 evolution is dominated by the g — > qq splitting (see section 4.1 and [86]). The compatibility 
of 6 R with zero for (3 < 0.3 thus implies that the ratio of the quark to the gluon density is 
similar in the diffractive and inclusive cases when considered at the same x values. Indeed, 
global fits to inclusive DIS data [41,78,79] yield gluon fractions of approximately 70% at low 
x, compatible with the results of section 5.4. At higher (3 > 0.3, where the DPDFs develop a 
more complicated structure (see figure 11) and the q — > qg splitting becomes important in the 
evolution (see section 4.1), b R becomes negative. 

The ratio defined in equation 18 can also be plotted as a function of x (or f3) with x F and 
Q 2 fixed. However, this results in a complicated dependence, which is driven by the high f3 
structure of the diffractive reduced cross section (figures 2-6a). 

The x dependence of the ratio of the diffractive to the inclusive cross section has been 
studied previously 10 at fixed M x [8,9, 11] rather than fixed x p . Using the diffractive data as 
presented in figure 8, the quantity 

is formed. Assuming proton vertex factorisation and neglecting contributions from sub-leading 
exchanges, the generalised optical theorem [85] predicts that this ratio is independent of Q 2 and 
depends only weakly on (3 — Q 2 /{Q 2 + M x ) and x ~ Q 2 /W 2 for sufficiently large M x . In 
models in which both the diffractive and the inclusive cross sections are governed by a universal 
pomeron [84], the remaining weak x dependence of the ratio arises due to the deviations from 
unity of the pomeron trajectory. 

The ratio defined in equation 20 is shown in figure 19 as a function of x in bins of fixed Q 2 , 
(3 and hence M x . In order to simplify the interpretation, data points are excluded if, according 
to the QCD fit in section 5, the sub-leading exchange contribution is larger than 10% or cr®^ 
differs by more than 10% from F®^ due to the influence of F^ 3 \ Only Q 2 and f3 values for 
which there are at least two remaining data points are shown. The ratio in equation 20 is indeed 
approximately constant throughout the full kinematic range, except at large j3 values (the low 
M x "non-triple-Regge" region in which [85] is not applicable). In particular, the x (and hence 
the W) dependence at fixed Q 2 , (3 and M x is strikingly flat, substantiating the conclusions 
of [8,9, 11]. 

Expressed in terms of Regge trajectories, the ratio of cross sections shown in figure 19 is 
proportional to x~ K , where k = 2(a F (t)) — a™ cl (0) — 1. Here ct p (t) is the effective pomeron 
trajectory for diffractive DIS and a™ cl (0) is the effective pomeron intercept governing inclusive 
scattering. Analysis of inclusive DIS data has shown that a™ cl (0) is not universal, but varies 
with Q 2 according to c^ cl (0) ~ 1 + 0.048 In (Q 2 / [0.292 GeV] 2 ) [87]. A prediction is overlaid 

10 The analyses in [8, 9, 1 1] differ from that described here in that the data were plotted as a function of W ~ 
y / Q 2 /x. In [8, 1 1], the results were also shown integrated over ranges in M x , whereas they are presented here at 
fixed M x values. 



25 



in figure 19, where it is assumed that the diffractive a F (0) depends on Q 2 in the same way as 
a mcl (0) and that a' = 0.06 GeV™ 2 (see section 5.2). Similar results are obtained with a' = 

IP V / IP v ' IP 

or a' p = 0.25. The normalisation of the prediction is obtained from separate fits to the data 
for each pair of f3 and Q 2 values. The description of the data is poor (x 2 = 876 based on 
statistical and uncorrelated systematic errors for 223 degrees of freedom). The ratio shown in 
figure 19 is also compared with a prediction where the inclusive data are described by the same 
Q 2 dependent a™ cl (0) of [87] and the diffractive x dependence is determined by the flux factor 
defined in equation 14, such that ct F (0) = 1.118 (equation 16) independently of Q 2 . A much 
improved description is obtained (x 2 = 254), with an x dependence of the ratio which changes 
slowly with Q 2 , being approximately flat in the region of Q 2 ~ 15 GeV 2 . 

The ratio of cross sections shown in figure 19 is incompatible with a simple Regge approach 
to 7*p scattering, where both the diffractive and the inclusive cross sections are driven by the 
exchange of the same pomeron trajectory, even if that trajectory is allowed to change with Q 2 . 
There is thus no simple relationship within a Regge model between the diffractive and inclusive 
cross sections. The flatness of the ratio of cross sections is natural if rapidity gap formation is a 
random process [14, 15] and has also been interpreted in colour dipole models [18]. 



8 Summary 

The reduced semi-inclusive diffractive cross section a^ 3 \x p , (3, Q 2 ) is measured for the pro- 
cess ep — > eXY under the conditions M Y < 1.6 GeV and \t\ < 1 GeV 2 for various fixed values 
of x p in the range 0.0003 < x p < 0.03. The data span nearly three orders of magnitude in Q 2 
from 3.5 GeV 2 to 1600 GeV 2 and cover the range 0.0017 < (3 < 0.8. In the best measured 
region, the statistical and point-to-point systematic uncertainties are at the level of 5% each, 
with an additional normalisation uncertainty of 6%. The kinematic dependences of the ratio of 
the diffractive to the inclusive reduced cross section are also studied. 

Up to small deviations at large x F and low (3, which are consistent with expectations from 
the presence of a sub-leading exchange, the f3 and Q 2 dependences of the diffractive data change 
only in normalisation at different x F values. This remarkable feature is compatible with a 
factorisable proton vertex. The variation of with x p can be expressed in terms of an 
effective pomeron trajectory with intercept 0:^(0) = 1.118 ± 0.008 (exp.) ig^io (model) if 
a' F = 0.06 Ip og GeV™ 2 is taken from HI Forward Proton Spectrometer data. The x (or x p ) 
dependence of the diffractive cross section at fixed f3 and Q 2 is similar to that of the inclusive 
cross section. The diffractive and inclusive x dependences cannot be interpreted with a single 
Q 2 dependent effective pomeron trajectory. 

The (3 and Q 2 dependences of ov^ 3 -* are interpreted in terms of diffractive parton distribution 
functions (DPDFs), obtained through an NLO DGLAP QCD fit. The DPDFs correspond to 
integrals over the measured M Y and t ranges and are valid in the region Q 2 > 8.5 GeV 2 , 
(3 < 0.8 and M x > 2 GeV. At low (3 < 0.3, the diffractive data exhibit a rather fast rise 
with increasing Q 2 at fixed x p and x. This rise of the diffractive cross section is very similar 
to that of the inclusive cross section at the same x values, implying that the ratio of quarks to 
gluons is similar in the diffractive and inclusive cases. The low f3 data give good constraints on 
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the diffractive quark singlet and gluon densities at low momentum fractions z, with combined 
experimental and theoretical uncertainties of typically 5—10% and 15%, respectively. The gluon 
density dominates the DPDFs in this region. At larger z, the diffractive quark density remains 
well constrained by the fit, whereas the sensitivity to the gluon density becomes increasingly 
poor. 

At high Q 2 , charged current scattering is used to test the assumptions on the quark flavour 
decomposition of the DPDFs. Total and single differential diffractive e + p charged current 
cross sections are measured and are well described by predictions based on the DPDFs ex- 
tracted from the neutral current data, though the current level of experimental precision (35% 
for the total cross section) is low. The DPDFs will provide important input to future tests of 
the factorisation properties of diffraction and to the prediction of cross sections for diffractive 
processes at HERA, the LHC and elsewhere. 
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Figure 2: The (3 and Q 2 dependences of the diffractive reduced cross section, multiplied by x p , at x p = 0.0003. In (b) the data are 
multiplied by a further factor of 3 l for visibility, with i as indicated. The inner and outer error bars on the data points represent the statistical 
and total uncertainties, respectively. Normalisation uncertainties are not shown. The data are compared with the reduced cross section at 
E p = 820 GeV derived from the results of 'HI 2006 DPDF Fit A', which is shown as a shaded error band (experimental uncertainties only) 
in kinematic regions which are included in the fit and as a pair of dashed lines in regions which are excluded from the fit. 
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Figure 3: The (3 and Q 2 dependences of the diffractive reduced cross section, multiplied by x F , at x p = 0.001. See the caption of figure 2 
for further details. 
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Figure 4: The j3 and Q 2 dependences of the diffractive reduced cross section, multiplied by x F , at x F = 0.003. In (a), the quantity 



y 2 /Y + ■ F L is also shown, as extracted from the 'HI 2006 DPDF Fit A'. Adding this quantity to the reduced cross section yields 
See the caption of figure 2 for further details. 
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Figure 5: The (3 and Q 2 dependences of the diffractive reduced cross section, multiplied by x p , at x p = 0.01. In (a), the contribution of the 
sub-leading exchange alone according to the 'HI 2006 DPDF Fit A' is also shown. The data with Q 2 < 90 GeV 2 (Q 2 > 200 GeV 2 ) were 
obtained with E p = 820 GeV (E p = 920 GeV). The fit results are shown for E p = 820 GeV. See the caption of figure 2 for further details. 
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Figure 7: (a) Measurements of the logarithmic Q 2 derivative of the reduced diffractive cross 
section at different fixed values of x p and (3, obtained by fitting equation 10 to the data. The 
derivatives are divided by the diffractive flux factor as defined in equation 14. (b) The loga- 
rithmic Q 2 derivative at x F = 0.01, divided by the diffractive flux factor and compared with 
the prediction of the 'HI 2006 DPDF Fit A'. The prediction is also decomposed into contribu- 
tions to the evolution from the splittings g — > qq ('Gluon driven evolution') and q — >• qg ('Quark 
driven evolution'). The inner and outer error bars represent the statistical and total uncertainties, 
respectively. Normalisation uncertainties are not shown. 
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and total uncertainties, respectively. Normalisation uncertainties are not shown. The data are 
compared with the results of the 'HI 2006 DPDF Fit A' for E p = 820 GeV, which is shown as 
a shaded error band (experimental uncertainties only) in kinematic regions which are included 
in the fit and as a pair of dashed lines in regions which are excluded from the fit. 
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Figure 9: Comparisons between subsets of the present data and results obtained by the direct 
measurement of the final state proton using (a) the HI Forward Proton Spectrometer (FPS) [6] 
and (b) the ZEUS Leading Proton Spectrometer (LPS) [10]. The FPS and LPS data are shifted 
to the Q 2 and (3 values shown using small translation factors and are multiplied by a further 
universal factor of 1.23 such that they correspond to M Y < 1.6 GeV. The inner error bars 
represent the statistical uncertainties and the outer error bars represent the combined statistical 
and systematic uncertainties. Normalisation uncertainties are not shown. 
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Figure 10: Comparison between a subset of the present data and the results of the ZEUS col- 
laboration, obtained by decomposition of the inclusive In M 2 distribution [11]. The ZEUS data 
are shifted to the Q 2 and j3 values shown using small translation factors and have been multi- 
plied by a further universal factor of 0.86 so that they correspond to M Y < 1.6 GeV. The inner 
error bars represent the statistical uncertainties and the outer error bars represent the combined 
statistical and systematic uncertainties. Normalisation uncertainties are not shown. 
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Figure 11: The total quark singlet and gluon distributions obtained from the NLO QCD 'HI 
2006 DPDF Fit A', shown at four different values of Q 2 for the range 0.0043 < z < 0.8, 
corresponding approximately to that of the measurement. The light coloured central lines are 
surrounded by inner error bands corresponding to the experimental uncertainties and outer error 
bands corresponding to the experimental and theoretical uncertainties added in quadrature. 
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Figure 12: Comparison on a linear z scale between the total quark singlet and gluon distribu- 
tions obtained from the 'HI 2006 DPDF Fit A' and the 'HI 2006 DPDF Fit B'. These two fits 
differ in the parameterisation chosen for the gluon density at the starting scale for QCD evo- 
lution. The DPDFs are shown at four different values of Q 2 for the range 0.0043 < z < 0.8, 
corresponding approximately to that of the measurement. For 'Fit A', the central result is shown 
as a light coloured central line, which is surrounded by inner error bands corresponding to the 
experimental uncertainties and outer error bands corresponding to the experimental and theo- 
retical uncertainties added in quadrature. For 'Fit B', only the total uncertainty is shown. 
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Figure 13: The Q 2 dependence of the fraction of the longitudinal momentum of the diffractive 
exchange which is carried by gluons according to the 'HI 2006 DPDF Fit A' , integrated over the 
range 0.0043 < z < 0.8, corresponding approximately to that of the measurement. The light 
coloured central line is surrounded by an inner error band corresponding to the experimental 
uncertainty and outer error band corresponding to the experimental and theoretical uncertainties 
added in quadrature. The central result from the 'HI 2006 DPDF Fit B' is also indicated. 
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or P, with error bars corresponding to the full experimental uncertainties. The bands show the 
result and experimental uncertainty from the standard fit in which a single parameter is used for 
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Figure 15: The cross section for the diffractive process e + p — > v e XY at E p = 920 GeV, shown 
differentially in (a) x p , (b) (3 and (c) Q 2 , for Q 2 > 200 GeV 2 , y < 0.9 and x p < 0.05. The 
data are compared with the predictions of the 'HI 2006 DPDF Fit A' to the neutral current 
data, obtained using the RAPGAP Monte Carlo generator. The contribution to the predictions 
from the sub-leading exchange is also shown. The data points correspond to average values 
of the differential cross sections over the regions shown by the histograms. The inner and 
outer error bars on the data points represent the statistical and total uncertainties, respectively. 
Normalisation uncertainties are not shown. 
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Figure 16: The ratio of the diffractive to the inclusive reduced cross section, multiplied by x p and shown as a function of Q 2 for fixed x 
and fixed x p = 0.0003, 0.001 and 0.003. The data are multiplied by a further factor of 3* for visibility, with % as indicated. The inner and 
outer error bars represent the statistical and total uncertainties, respectively. Normalisation uncertainties are not shown. The results of fits of 
a linear dependence on logQ 2 (equation 19) to the data are also shown. 
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Figure 17: The ratio of the diffractive to the inclusive reduced cross section, multiplied by x F and shown as a function of Q 2 for fixed x and 
fixed x F = 0.01 and 0.03. See the caption of figure 16 for further details. 
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Figure 18: The logarithmic Q 2 derivative of the ratio of the reduced diffractive cross section 
to the reduced inclusive cross section at different fixed values of x p and (3 (see equation 19). 
The derivatives are divided by the diffractive flux factor as defined in equation 14. The inner 
and outer error bars represent the statistical and total uncertainties, respectively. Normalisation 
uncertainties are not shown. 
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Figure 19: The ratio of the diffractive to the inclusive reduced cross section, multiplied by 
(1 — (3) x p . Data at /3 = 0.9 are multiplied by a further factor of 5 for visibility. Data which 
are significantly influenced by the sub-leading exchange or the longitudinal structure function 
according to the NLO QCD fit are excluded. The remaining data are compared with models 
motivated by Regge phenomenology as described in the text, for which the diffractive and 
inclusive effective pomeron trajectories are the same ('a p = a inc1 ') and for which they are 
different ('q p ^ a> c1 '). 
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Table 4: Results for x p o r at fixed Q 2 , (3 and x p (columns 1-4) using data with SPACAL electrons and E p = 820 GeV. Columns 5-7 
contain the percentage statistical, systematic and total uncertainties. The remaining columns contain the contributions to the systematic 
uncertainty from sources which are uncorrelated between data points (5 unc ) and the 10 correlated sources leading to the largest uncertainties. 
These are the LAr hadronic energy scale (Si ar ), the SPACAL electromagnetic energy scale (5 e i e ), the scattered electron angle measurement 
(Sg), the calorimeter noise treatment (5 no i se ), reweighting the simulation in x F (5 Xp ) and (3 (dp), the background subtraction using DJANGO 
(8b g ), the plug energy scale (5pi ug ), reweighting the simulation in Q 2 (5q2) and the SPACAL hadronic energy scale (5 spa ). Minus signs appear 
for these systematics if the shift in a variable is anti-correlated rather than correlated with the shift in the cross section. The table continues 
on the next 7 pages. 
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-2.3 


1.7 


0.6 


0.03 


200.0 


0.27 


0.0283 


5.3 


12.3 


13.4 




0.1 


0.5 


-0.8 


0.8 


1.7 


-0.9 


3.8 


-2.1 


1.4 


0.2 


0.03 


200.0 


0.43 


0.0309 


5.5 


12.2 


13.4 




-0.5 


1.2 


-0.7 


0.1 


1.5 


-1.0 


3.7 


-2.1 


0.9 


-0.1 


0.03 


200.0 


0.67 


0.0297 


7.2 


13.4 


15.2 




-1.0 


2.8 


-0.9 


-1.2 


0.8 


-0.6 


5.8 


-2.6 


0.2 


0.5 


0.03 


400.0 


0.27 


0.0322 


7.6 


13.4 


15.4 




-0.4 


0.6 


-0.2 


0.1 


1.2 


-0.8 


6.9 


-1.8 


1.3 


-0.3 


0.03 


400.0 


0.43 


0.0293 


6.9 


12.2 


14.0 




-0.7 


-0.4 


-0.1 


-0.3 


0.9 


-0.8 


4.1 


-2.2 


0.7 


-0.2 


0.03 


400.0 


0.67 


0.0289 


8.0 


13.5 


15.6 




-0.9 


1.2 


-0.7 


-0.9 


1.1 


-0.9 


6.6 


-2.1 


0.5 


-0.4 


0.03 


800.0 


0.43 


0.0391 


13.1 


13.7 


19.0 




-0.9 


2.4 


-1.2 


-0.3 


0.1 


-0.5 


7.1 


-2.1 


0.2 


-0.4 


0.03 


800.0 


0.67 


0.0228 


14.6 


14.6 


20.6 




-0.1 


-0.9 


0.3 


-0.5 


0.2 


-0.5 


9.1 


-2.1 


0.1 


-0.2 


0.03 


1600.0 


0.67 


0.0214 


27.9 


15.5 


31.9 




-1.6 


-5.7 


-3.2 


-1.7 


-0.4 


-0.1 


7.5 


-2.6 


-0.2 


-0.3 



Table 5: Results for x p a® at fixed Q 2 , (3 and x p (columns 1-4) using data with LAr electrons and E p = 920 GeV. Columns 5-7 contain 
the percentage statistical, systematic and total uncertainties. The remaining columns contain the contributions to the systematic uncertainty 
from sources which are uncorrelated between data points (5 unc ) and the 10 correlated sources leading to the largest uncertainties. These 
are the LAr hadronic energy scale (Si ar ), the LAr electromagnetic energy scale (S e i e ), the scattered electron angle measurement (So), the 
calorimeter noise treatment (5 noise ), reweighting the simulation in x p (8 X ) and (3 (6/3), the background subtraction using the non-diffractive 
RAPGAP simulation (5b g ), the plug energy scale (Spi ug ), the 1 — (3 reweighting of the simulation (Spm) and the contribution to the hadronic 
energy from charged particle tracks (5 tra )- Minus signs appear for these systematics if the shift in a variable is anti-correlated rather than 
correlated with the shift in the cross section. 



x p range 


(x p ) 


da/dx p [pb] 
Q 2 > 200 GeV 2 , y < 0.9, x p < 0.05 


Vstat 
(%) 


(%) 


(%) 


0.005 - 0.016 
0.016 - 0.05 


0.01 
0.03 


15 
7.7 


59 
37 


25 
18 


64 
41 



g 2 range 

[GeV 2 ] 


(Q 2 ) 

[GeV 2 ] 


da/dQ 2 [pb/GeV 2 ] 
Q 2 > 200 GeV 2 , y < 0.9, x p < 0.05 


Ostat 
(%) 


Osys 
(%) 


(%) 


200 - 560 
560 - 2240 


350 
1150 


6.0 • 10" 4 
1.0 • 10" 4 


50 
40 


21 
18 


54 
44 



/3 range 


<0> 


da/dp [pb] 
g 2 > 200 GeV 2 , y < 0.9, x p < 0.05 


O s tat 
(%) 


Osys 
(%) 


(%) 


0.3 - 0.5 


0.40 


0.58 


62 


27 


68 


0.5-0.8 


0.65 


0.44 


44 


21 


48 


0.8-1.0 


0.90 


0.17 


79 


17 


81 



Table 6: Measurements of the cross section for the diffractive charged current process e + p — > 
z/ e XV for g 2 > 200 GeV 2 , y < 0.9 and < 0.05 with £ p = 920 GeV, differential in x p , Q 2 
and p. The differential cross sections correspond to average values over the ranges shown. The 
percentage statistical, systematic and total uncertainties are also given. 
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